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The spin- valve complex magnetoimpedance of symmetric ferromagnet /normal metal/ferromagnet 
junctions is investigated within the drift-diffusion (standard) model of spin injection. The ac 
magnetoresistance — the real part difference of the impedances of the parallel and antiparallel mag- 
netization configurations — exhibits an overall damped oscillatory behavior, as an interplay of the 
diffusion and spin relaxation times. In wide junctions the ac magnetoresistance oscillates between 
positive and negative values, reflecting resonant amplification and depletion of the spin accumula- 
tion, while the line shape for thin tunnel junctions is predicted to be purely Lorentzian. The ac 
spin-valve effect could be a technique to extract spin transport and spin relaxation parameters in 
the absence of a magnetic field and for a fixed sample size. 

PACS numbers: 72.25.Ba, 72.25.Rb, 85.75.-d 



Electrical spin injection from a ferromagnetic (F) to 
a nonmagnetic {N) conductor is essential for spintron- 
ics [H, Ql- Predicted byAronov Q and first realized by 
Johnson and Silsbee, [^-|6[, it is now a well established 
concept. A biased ferromagnetic/nonmagnetic junction 
generates a noncquilibrium spin accumulation within the 
spin diffusion length at the interface, building a noncqui- 
librium resistance 0, In an FNF junction this 
noncquilibrium resistance gives rise to the difference in 
the junction resistances for parallel (P) and antiparallel 
(AP) magnetization orientations of the F regions — the 
giant magnetoresistance (GMR) 0, [l^l- Drift-diffusion 
theory along with the spin accumulation concept success- 
fully describes magnetoresistance effects in charge neu- 
tral [ill [l3 ^^^^ space-charge systems (l3. 14|. 
enabling one to obtain relevant spin-related materials pa- 
rameters 15[, such as the spin relaxation times. 

Recently, Rashba has generalized the spin-polarized 
drift-diffusion theory to the alternating current (ac) 
regime [1, [l^. We apply this theory and investigate 
the complex impedance Z{u!) of symmetric FNF junc- 
tions. Wc show that the real part of the spin- valve mag- 
netoimpedance (we call it here ac magnetoresistance) 
AZ ~ Zap — Zp of the junctions exhibits damped oscil- 
lations as a function of frequency. The oscillation period 
is given by the diffusion time through the normal layer. 
In mesoscopic junctions (of sizes up to the spin relaxation 
length Lg), the ac magnetoresistance can be negative at 
experimentally accessible frequencies, meaning that the 
antiparallel configuration has a lower ac resistance than 
the parallel one. The negative ac magnetoresistance is 
a consequence of a resonant spin accumulation effect, 
namely a resonant spin amplification in the P configu- 
ration and a resonant spin depletion in the AP one. In 
nanoscale junctions (with sizes much less than Lg), with 
tunnel contacts, the oscillation period is large, leaving a 
nice Lorentzian profile with the width of the spin relax- 
ation rate. A one-parameter fit to the line shape (either 
damped oscillator or Lorentzian) determines the spin re- 
laxation time T,. 



We present the ac spin-valve effect as an alternative 
to other methods that measure of nonmagnetic con- 
ductors, such as the conduction electron spin resonance, 
spin pumping, or the Hanle effect, which require mag- 
netic fields, or to the dc spin injection method (in verti- 
cal or lateral geometries), which requires studying various 
sample sizes (distances to electrodes) to extract the spin 
diffusion length [13. In a sense the ac spin- valve effect 
is similar to the Hanle effect, which is widely used to 
fin d sp in relaxation times in metals and semiconductors 
[3, Il8l |. but the role of the magnetic field is taken by the 
frequency; in the Hanle effect too the signal in general 
oscillates as a function of magnetic field, with a modified 
Lorentzian shape in the diffusive regime [2|. 
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FIG. 1: Scheme of an FNF spin valve. The spacer A'^ region 
has width d and the sizes of the ferromagnetic electrodes are 
assumed greater than the corresponding spin diffusion length 
Ls- In the dc regime the parallel configurations results in 
smaller spin accumulation (dashed line) than in the antipar- 
allel one, demonstrated by the positive dc spin-valve magne- 
toresistance. In the ac regime, this can be reversed (solid): at 
certain frequency ranges there can be a resonant spin ampli- 
fication in the parallel and spin depletion in the antiparallel 
configuration, resulting in a negative ac magnetoresistance. 



A microwave measurement of Tg in the absence of a 
magnetic field may be important as in many conduc- 
tors the spin relaxation time depends strongly on it; a 
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striking case is aluminum in which decreases by an 
order of magnitude as the magnetic field increases from 
0.05 to 1.3 T [l^. Still, Ts obtained by spin resonance 
tend to be, for a given temperature, much greater than 
that obtained from tr ansp ort techniques, as catalogued 
for Al and Cu in Ref. |20l The case of Au is even more 
striking, as spin resonance shows that at low tempera- 
tures the ratio of Tg to the momentum relaxation time is 
about one, while transport techniques predict the ratio 
to be about 100 2lJ, [221; at room temperature, at which 
phonons are relevant, the ratio is about 10, as measured 
by spin pumping (23j which requires both magnetic field 
and nanoscale transparent junctions. For extracting bulk 
spin relaxation times it may be preferable to work with 
tunnel contacts and mesoscopic samples, so that spin re- 
laxation is not strongly influenced by the interface and 
surface effects. (Various techniques for measuring as 
well as useful data are given in the review Ref. Il5|.) The 
ac spin-valve method could potentially explore nano and 
mesoscopic spin valves, in both vertical and lateral ge- 
ometries, at no magnetic field applied to the normal con- 
ductor, and provide the spin relaxation times at a fixed 
sample size [27[. 

We consider a symmetric FNF junction as compris- 
ing two FN junctions in series, see Fig. [TJ Each FN 
junction has a contact (c) region with a spin-dependent 
conductance; otherwise spin is assumed to be preserved 
at the contact. The spin valve dc magnetoresistance 
AR. = AZ{uj = 0) of a symmetric FNF junction, whose 
N region has width d and the F regions have widths much 
greater than the spin diffusion lengths, can be expressed 
analytically within the drift-diffusive regime [l|,[l,|lB: 24|. 
This dc formula has a straightforward extension to the 
harmonic ac regime, and we write the complex magne- 
toimpedance as 



AZ{uj,d)^ (1) 

8rjv(cj)[rF(c^) P,^ + r, PsJ' e''/^-"''") 

[rpiuj] +rc + riv(a;)]2 e2d/L,„(tu) _ [rp{uj) + Vc - rAr(w)]2 

by indicating the complex frequency-dependent quanti- 
ties (labeled by the region and F), 



We analyze the spin- valve impedance, based on 
Eq. ([T]), for a realistic model Py/Cu/Py junction, with 



Ls{uj) = Lg I Vl - iWTs, 

r(uj) = r I \/\ — iLOTg. 



(2) 
(3) 
(4) 



Here is the spin relaxation time, Lg = VDtI is the 
spin diffusion length, D is the diffusivity, and r — Lg/cr is 
the effective resistance, with a denoting the conductivity. 
The effective contact resistance is Vc = {T,^ + T,i)/4'Ef'Ei, 
with the contact conductance of spin A. Finally, Pa-^^ 
and Pe^ denote the spin polarization of the conductivity 
and conductance of the F and contact region, respec- 
tively. Driving ac is assumed to be harmonic with the 
angular frequency lo — 27r/, i.e. j{t) oc e~*'^*. 



the following experimentally obtained data 
the temperature T = 4.2 K: Lgpf = 1 /im, r^jv 
Dn = 238cm2s-\ rjv 
TsF = 0.6 ps. Dp 

Pa. 
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at 
42 ps, 

= 14ffim^, LsF = 5.5 nm, 
rp = 0.42 fOm^, 
0.22. For the contact characteristics we employ 



[25|: rc{— rp) = 0.5ff2m^ and Ps^ = 0.4, so the con- 
tact interface is generic, neither tunnel nor transparent. 
The specific spin resistivities rp, and Tc and hence 
the spin valve AZ{uj, d) are evaluated for a unit cross 
section. In the experiment one divides these resistivities 
by the actual conductor cross sections, which could be 
10-3- l^m^. 

Figure [2] presents the calculated magnetoresistance. In 
Fig.HKb) we show the dc magnetoresistance as a function 
of the d. With increasing d the magnetoresistance expo- 
nentially decreases, as the injected spin accumulation is 
damped. The plot in Fig. ^c) shows the ratio of the 
ac to the dc magnetoresistance, Re[AZ{f,d)]/ATZ{d), as 
a function of d and frequency / = a;/27r. For a given 
d, the ac magnetoresistance oscillates as a function of /, 
between positive and negative values. In Fig.[2][a) the os- 
cillations are shown for d = 4 /im. The negative peaks are 
considerable fractions (tens of percents) of the dc values. 
On the d-f -plot the oscillations show hyperbolic stripes. 
For thin samples, the dependence on / is rather weak for 
this generic junction. We will sec below that for tunnel 
junctions the dependence becomes Lorentzian. 
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FIG. 2: (color online) ac spin-valve effect in a model 
Py/Cu/Py junction, (a) Calculated ac/dc ratio of the spin- 
valve magnetoresistance as a function of / = u}/2n for d = 
4 fj.m. (b) Calculated dc spin- valve magnetoresistance as a 
function of d. (c) Calculated ac/dc ratio of the spin- valve 
magnetoresistance as a function of d and driving frequency /. 
The visible light and dark bands of equal signs are separated 
by the node lines, Ile[AZ{f, d)] — 0. 



To be specific, consider d 



sN 



1 fiia. The ac 
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magnetoresistance remains positive for / < 34.8 GHz. 
Further increase in the driving frequency leads to a neg- 
ative ac spin-valve magnetoresistance: Re(ZiZ) < 0. For 
d = 3LsN = 3 /im the spin- valve magnetoresistance re- 
mains positive up to the frequency / ~ 6 GHz, then it 
becomes negative for 6 GHz < / < 26.1GHz. There 
should be more oscillations observable at larger values of 
d, but at the cost of exponentially reducing the magni- 
tude, see Fig. [2lja). We will see below that the relevant 
time scale parameter for the oscillations is the diffusion 
time through the spacer layer. For our model junction, 
a reasonable parameter range for measuring the ac oscil- 
lations would be the sample sizes Lsn ^ d ^ '^Lsn- The 
involved frequency, / = a;/27r, ranges are 1 GHz— 50 GHz, 
experimentally well accessible. 

Mathematically, the spin valve oscillations appear nat- 
urally. The real dc transport parameters become in the 
ac case complex, see Eqs. ©-([I]). The imaginary part 
of d/LsN{^) gives rise to the ac exponential Qd./LsN[<^) 
in Eq. ([T]) with the trigonometric character and hence 
a certain oscillatory behavior of the complex spin valve 
impedance AZ{u). For the frequencies uj <^ T^Ni"^ '''If) 
the imaginary part oi d/ Lsn{^) plays no role, see Eq. ([3|). 
The ac magnetoresistance exhibits changes on the scales 
of the relaxation rate 1 /ts or the diffusion rate through 
the spacer. These provide the practical limit for the use 
of microwaves in the experiment. 

We now give a qualitative picture of the predicted os- 
cillatory behavior, including the negative ac spin-valve 
magnetoresistance. First, we show that the spin valve 
impedance AZ{uj) is related to the contact values of the 
spin accumulations in TV, for P and AP configurations. 
From the standard spin injection model for a symmetric 
FNF junction we derive the following formula [26| : 



Jit) 



AZ{uj) . 



(5) 

Here fJ'^lf^^ (c, t) represent the actual nonequilibrium spin 
accumulation in the N spacer for P and AP configura- 
tions respectively, at the left FN contact interface c (see 
Fig. [ij and j(t) is the driving harmonic ac. To under- 
stand the ac magnetoresistance oscillations, one needs to 
look at the contact spin accumulation only. 

The qualitative picture is in Fig. [31 which shows P and 
AP configurations at three times t = 0, t = Tn/A, and 
t = Tn/2, where Tn = Tsn d'^ / LIj^ = d^/DN is the dif- 
fusion time through N . The resonant spin amplification 
and depletion effect happens if the driving current j{t) 
is close to the N spacer diffusion time T/v, this case is 
shown in Fig. |3l 

At time t = the current j is negative and electrons 
are injected from the left and extracted to the right elec- 
trodes, leaving behind positive and negative spin accu- 
mulations, indicated in Fig. |3] by diffusive packets (the 
sample is locally charge neutral, only spin is redistributed 
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FIG. 3: (color online) Mechanism for the resonant amplifica- 
tion and depletion of the spin accumulation in an ac-driven 
FNF junction for parallel (P) and antiparallel (AP) configu- 
rations. The solid lines within the horizontal of the N spacer 
represent actual profiles of the nonequilibrium spin accumula- 
tion fj^lf^^ {x,t), which correspond to the harmonic ac signal 
j{t) shown on the left. The injected and extracted spin pack- 
ets and their diffused and spread positions are shown at three 
distinct times t = 0, t = Tm/A, and t = Tn/2, where Tn 
is the characteristic diffusion time across the N conductor. 
Frequency of the driving ac is close to 1/Tn and tones of the 
packets correspond to F conductors which initially emitted 
them. 



nonuniformly) . The dynamics of these spin packets is 
governed by diffusion and relaxation, but by not bias 
voltage. This is because in the N spacer, there is no 
spin-charge coupling and spin and charge transports are 
decoupled, see 13- At time t = T/v/4 the current van- 
ishes, j = 0, as well as the spin injection and extraction. 
In the meantime the spin packets diffusively spread and 
reach the center of the N spacer. At t = Tn/2 the spin 
packets reach the other contacts. Now the current is fully 
reversed: in the P configuration the new spin packet is 
injected at the right electrode, amplifying the initial in- 
jected spin packet that has traveled from the left. In the 
AP configuration, the new spin packet of the opposite 
sign is injected at the right, depleting the initial injected 
spin. Similarly at the left electrode. The left contact dif- 
ference lJ>^ff{c) — M^7v('^) Sit t ~ Tn/2 becomes negative, 
the actual current j > and, according to Eq. ([5]), we 
get negative ac magnetoresistance, Rc{AZ) < 0. 

The resonance condition, ujTn ~ tt, is equivalent 
to LsN/d « y/uTsN/TT. In practice, to see the nega- 
tive ac magnetoresistance one prefers d ~ Lsn, so that 
ujTsN ~ Ti", which is the microwave regime. If d > Lsn, 
as in our model shown in Fig. [21 then the oscillations can 
be observed at lower frequencies, but at the cost of de- 
creasing the magnitude of the ac magnetoresistance due 
to spin relaxation. This need not be an issue with tunnel 
contacts, as the precision of measuring higher resistances 
is higher. On the other hand, no oscillations (within the 
GHz regime) should be seen for nanoscale junctions, for 
d <^ LsN- We will show below that this important regime 
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gives a Lorentzian profile. 

We now turn to the case of a junction with tunnel 
contacts, rp,rN <ti r^- (In general, using tunnel barri- 
ers allows to adopt the junction resistance and size in- 
dependently, providing maximal flexibility in the device 
design.) In this case Eq. ([T]) reduces to 



AZ 



4 TAT 



Pi 



(6) 



where Lgpf = \^DnTsn- A single parameter {tsn, know- 
ing d and Dn) fit of a measurement of the w dependence 
of the tunnel spin valve impedance (relative to the dc 
value) to Eq. ([5]) can determine the spin relaxation time 
of the normal region. The shape is illustrated in Fig. 3] 
Since Tsf is typically an order or two magnitudes smaller 
than TsNi the ac effects do not play significant in the F 
electrodes. 

For d > |Lsjv(w)| we can approximate AZ as follows: 




FIG. 4: Calculated ac/dc ratio of tlie FNF magnetoresis- 
tances for three different widths d of the A'" spacer with tun- 
nel contacts to F. The solid line, d — 0.1 /im, represents 
a Lorentzian line shape with its half-width determining the 
spin relaxation time Tsn- The dotted line, d = 2 /^m, shows 
the first oscillation, with negative ac magnetoresistance. The 
dashed line, d = 1 fim, is the intermediate case. The values 
for the parameters are as in Fig. [21 but with a greater tunnel 
resistance, Tc = 5nQm^. 



AZ{uj, d, T. 



sN) 



-d 



DnTsN 



lUJTsN 



(7) 



Suppose we know the experimental value of the frequency 
Wo at which Iic[AZ{ujQ, d)] vanishes. As an alternative 
to the fitting, the spin relaxation time can be given by 
the equation: 



1-h vr+^^ 



Q'sN 



tan 



dojQ 



N 



TsN 



1 + Vl + wgr^ 



(8) 



N 



which can be solved for Tsn with simple numerics. 

In the opposite important case of ci <C isAr, Eq. 
becomes a Lorentzian: 



AZ{tu,d,TsN) ~ 4 



I'N P%^ 1 + iuTsN 



(9) 



The half-width frequency at which 

Re[AZ{uJif2,d)] = ^Z\72.(c?) determines the spin 
relaxation time according to Ts^v = l/wi/2. This 
Lorentzian shape is rather robust for tunnel junctions, 
illustrated in Fig |4l which also shows an intermediate 
case of d w PsN ■ 

For transparent contacts [vc <C rp^rj^) and nanoscale 
junctions, d <^ L^n^ the magnetoimpedance is 
AZ{u:,d) « 2rp{oj)P^^, the square root of the 
Lorentzian, but with the width of I/t^f- The shape is 
therefore featureless in the microwave regime, unless the 
ferromagnetic contacts have a relatively large spin relax- 
ation time, in which case t^f could be determined from 
the measured shape. In mesoscopic junctions oscillations 
should be visible. 

In summary, we have presented a simple but robust 
theory of the ac spin-valve effect in symmetric FNF 
junctions and predict negative ac magnetoresistance due 



to resonant amplification and depletion of the spin accu- 
mulation in the normal metal region. The oscillating line 
shape allows a single-parameter (spin relaxation time) 
fitting for mesoscopic and nanoscale spin valves; for the 
latter a Lorentzian shape should be seen with tunnel con- 
tacts. 
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